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ABSTRACT 


A continuous-resistance,  continuous- 
capacitance  analog  model  has  been  constructed 
and  used  to  simulate  a fire-exposed  T-beam 
structure.  The  model  was  arranged  to  permit 
study  of  both  one  and  two-dimensional  transient 
heat  flow  and  to  allow  for  Newtonian  heat 
losses  to  the  cool  surroundings.  Comparison 
of  the  results  with  those  of  a numerical  solu- 
tion of  the  problem  show  agreement  within  11$, 
based  on  the  value  of  the  input  signal.  Con- 
sidering the  size  limitations  of  the  electrical 
model,  edge  effects,  and  the  variation  in  its 
electrical  properties,  the  agreement  is  con- 
sidered reasonable. 

1.  Introduction 

The  theory,  construction  details  and  performance  of  an 
electronic  instrument  which  permits  solution  of  transient  heat 
flow  problems  by  direct  electrical  analogy  has  been  previously 
described  £l7  and  the  instrument  has  been  used  to  solve  a num- 
ber of  one-dimensional  problems.  Since  the  method  involves 
the  use  of  "lumped"  resistance  and  capacitance  elements,  the 
extension  to  two  (or  three)  dimensions  could  be  accomplished  by 
the  use  of  discrete  (i.e.  lumped)  circuit  elements  arranged  in 
a two-  (or  three-)  dimensional  array.  Because  of  simplicity, 
space  considerations  and  electronic  circuitry  problems  associ- 
ated with  "fast-time"  instruments  of  this  type,  the  use  of  dis- 
tributed parameter  elements  would  present  distinct  advantages. 
Furthermore,  this  would  eliminate  the  "lumping"  error  resulting 
from  the  use  of  a limited  number  of  discrete  elements;  however, 
other  errors  associated  with  non-uniformity  of  the  conductive 
medium  would  generally  be  introduced. 

Conductive-sheet,  conductive-solid  and  conductive-liquid 
analogs  find  considerable  use  in  simulating  two-and  three- 


dimensional  fields  governed  by  Laplace's  equation,,  eo  go  steady- 
state  heat  flow  For  transient  heat  flow  involving  poten- 

tial energy  storing  elements,  however,  previous  use  appears  to 
have  been  limited  primarily  to  continuous-resistance,  lumped- 
capacitance  analog So  Apparently,  the  first  published  descrip- 
tion of  a method  utilizing  a distributed  resistance  and  capac- 
itance analog  was  described  by  Fatt  £~37° 

This  report  presents  data  obtained  by  use  of  a continuous- 
resistance,  continuous-capacitance  analog  simulating  a fire- 
exposed  T-beam  structure  and  compares  the  results  with  corres- 
ponding solutions  obtained  by  means  of  finite  difference  compu- 
tations Z*  47;  it  is  intended  to  demonstrate  the  feasibility  of 
using  such  a model  for  the  approximate  solution  of  one-  and 
two-dimensional  transient  heat  flow  problems » 

2o  Analog  Model  System 

An  ideal  two-dimensional  analog  consists  of  a layer  of 
dielectric  material  of  uniform  capacity,  to  one  surface  of  which 
is  applied  a layer  of  conductive  material  of  uniform  resistivity, 
the  other  surface  being  a perfectly  conducting  layer  connected 
to  grounds  Ideally  both  the  dielectric  and  the  resistance  layer 
should  exhibit  the  same  degree  of  homogeneity  and  isotropicity  as 
the  thermal  body  being  simulated 0 Of  the  more  readily  available 
dielectric  materials,  barium  titanate  was  chosen  even  though  it 
is  known  to  be  nonisotropic  to  some  extent . A 4 by  4 by  1/8  in0 
barium  titanate  sheet  containing  a painted  resistance  layer  on 
one  surface  and  silver  foil  on  the  other  surface  was  specially 
fabricated  by  the  Centralab  Division  of  Globe  Union  Mfg0  Go0  in 
February  1958o  This  was  the  largest  size  piece  which  could  be 
conveniently  produced  at  the  time;  however,  the  thickness  of  the 
dielectric  was  maintained  to  a tolerance  of  + ChOOl  in0 

Direct  measurements  were  made  of  the  resistance  at  several 
positions  across  the  entire  surface  and  these  are  listed,  in  Table 
1.  The  non-uniformity  of  the  resistance  layer  at  the  measuring 
point  locations  was  estimated,  by  means  of  voltage  measurements, 
to  be  about  ± 10  percents  The  time  constant  (RC)  measurements 
listed  in  Table  1 were  obtained  using  the  circuit  shown  in  Figo  1, 
A time  constant  of  7*+°7  msec,  is  well  suited  for  measurements  by 
conventional  cathode  ray  oseilloscopeso  From  these  measurements, 
the  capacitance  is  0o000h6h  mfd/cnr  which  corresponds  to  a 
dielectric  constant,  neglecting  edge  effects,  of  approximately 
l660o 
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Table  1.  Measured  Resistance  and  RC  Product 


R 

RC 

X-direction 

meg 

sec 

1 in. 

1.60 

.0752 

2 in. 

1.40 

3 in. 

1.40 

.074-7 

4 in. 

1.80 

.0770 

Y-direction 

1 in. 

2 in. 

1.60 

.0760 

3 in. 

.0741 

4 in. 

.0712 

Averages 

1. 56 

.0747 

The  analog  model  was  arranged  so  that  probes  could  be 
positioned  on  its  conductive  surface  through  vertically 
aligned  hole's  in  two  sheets  of  clear  plastic  placed  as  shown 
in  Fig.  2.  Although  a measuring  grid  of  this  type  was  selected 
for  convenience,  the  system  is  basically  continuous  and  does  not 
restrict  measurements  to  any  finite  number  of  space  points.,  Be- 
cause of  symmetry,  one-half  of  the  T-beam  may  be  represented  by 
the  inverted  L-section  of  Figure  3°  A series  of  19  "fixed" 
probes  were  placed  at  selected  points  and  were  spring-loaded  to 
ensure  good  contact  with  the  surface . The  tips  of  the  0.075  in. 
dia.  probes  were  beveled  to  provide  a contact  area  of  approxi- 
mately 0.03  in.  dia.  Each  of  the  probes  was  connected  to  an  in- 
dividual discharge  circuit  and,  through  a manual  selector  switch,  . 
to  the  high  impedance  measuring  circuit.  One  additional  "movable" 
probe  was  used  to  make  measurements  at  any  other  grid  point  de- 
sired. The  line  representing  the  fire-exposed  surface  was 


outlined  on  the  resistive  paint  surface  using  highly  conductive 
silver  paint  and  then  covered  with  a sheet  of  aluminum  foil 
formed  to  the  same  outline » The  lower  plastic  sheet  was  fas- 
tened down  firmly  assuring  good  contact  at  the  aluminum  foil- 
silver  paint  interface o A lead  connected  to  the  aluminum  foil 
electrode  was  fed  a voltage-time  signal  derived  from  a photo- 
former-type of  signal  generator  which  simulated  the  shape  of  the 
prescribed  temperature-time  curve . Another  lead  from  the  silver 
foil  on  the  lower  surface  of  the  model  was  connected  to  the  in- 
strument ground . 


The  problem  chosen  for  study  was  one  for  which  numerical 
solutions  by  means  of  a digital  computer  were  available  £47 „ 

The  lower  surfaces  (ceiling)  of  a long  horizontal  slab  and  beam 
structure  (T-beam)  are  subjected  to  the  standard  fire  exposure 
defined  for  fire  endurance  tests  £57  while  the  upper  surface 
(floor)  loses  heat  by  radiation  and  Newtonian  cooling  to  an 
ambient  maintained  at  the  initial  temperature  of  the  structure . 
Although  consideration  was  given  in  ref*  4 to  a wide  range  of 
physical  and  thermal  parameters,  this  study  was  concerned  only 
with  the  following: 


3 ° Two-dimensional  T-beam  Problem 


3°1  Thermal 


Thermal  conductivity: 


Thermal  diffusivity: 


Concrete  density: 


Beam  width: 


Beam  to  slab  ratio: 


Slab  thickness: 


Specific  heat: 


4o2 

4 in*  (10ol6  cm) 

8 in*  (20*32  cm) 

140  lb/ft^(2<.24  g/cm^) 

0*577  Btu/lb  F (2®4i 
watt  sec/g  C) 

0*905  Btu/hr  ft  F 
(0*0156  watt/cm  C) 

0*0112  f t2/hr  (0o 00289 


cm^/sec ) 
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The  loss  of  heat  from  the  unexposed  upper  surface  by  radi- 
ation and  convection  was  assumed  to  take  place  according  to  the 
relation 

V3 


q-  = 0o  275  t 


+ 0.173 


rr**4  -i 

o 

-p  1 

o 
o 
1 — 1 

n 

o 
o 
1 — 1 

p 

where  q = rate  of  heat  loss,  Btu/hr  ft 

t = temperature  rise  above  initial,  deg  F 
tQ=  initial  temperature,  68  F 


3o2  Electrical 

Analogous  thermal  and  electrical  elements  for  both  the  one- 
dimensional and  two-dimensional  cases  are  listed  in  Table  2. 


Table  2.  Analogous  Thermal  and 

Electrica] 

L Element So 

One- Dimensional 
Case 

Two-Dimensional  Case 
(Conductive  Sheet 
Analog ) 

Element 

Symbol 

Unit 

Symbol 

Unit 

Thermal  Resistance 

cm2C 

i 1 

cmC 

k 

watt 

k 

watt 

Electrical  Resistance 

R 

megohm 

R 

megohm/ 

square 

Thermal  Capacitance 

psxt 

watt-sec 

psAt 

watt-sec 

cm^C 

cm  C 

Electrical  Capacitance 

C 

cAe 

p fd 

In  addition  to  the  capacity,  resistivity  and  time-scaling  factors 
in  the  one-dimensional  case,  a linear  scale  factor  X = L+./Le 
is  introduced  in  the  two-dimensional  case  to  account  for  the  use 
of  a small-scale  electrical  model o Based  on  the  dimensions  and 
the  measured  electrical  values  of  the  analog  model,  the  values 
of  the  scaling  factors  (in  Cog„s.  units)  were  as  follows: 
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m =.  1 = 1 - 4l.l 

kR  (.0156) (1.56) 

n = ps  ,Lt  x2  = (2.24) (2.41)  (5) 2 = 2.91  x 10^ 
c VL  ; (.000464) 

9 " 1 

x = Lt  = 4 = 5 

L 0.8 

e 


For  a three-hour  fire  exposure  duration,  each  time  division 
on  the  read-out  oscilloscope  was  chosen  to  represent  4 minutes 
of  thermal  time  and  the  oscillator  was  therefore  set  for  a 
frequency  of  k 

mn  cb  (4l . 1)  (2. 91  x 10p)(l)  ^ kCe 

f = tt  = 4 x 60  ' " ' 

To  simulate  convective  and  radiative  heat  losses,  it  is 
common  to  use  non-linear  resistive  elements  [1,6].  To  maintain 
the  continuous,  geometrical  nature  of  a two-dimensional  analog 
would  require  the  use  of  a distributed  voltage-sensitive  resis- 
tive film.  In  the  absence  of  such  a material,  it  becomes  neces- 
sary to  employ  discrete  elements,  so  that  a combined  geometric 
and  network  analog  results.  The  size  of  the  connector  contact 
is  important  since,  if  too  large,  there  is  the  tendency  to  form 
a short  circuit  between  adjacent  contacts,  whereas  if  too  small, 
the  flow  lines  tend  to  crowd  [71.  Because  of  a lack  of  non- 
linear elements  of  the  desired  properties  and  since  simple  re- 
sistors may  approximate  non-linear  resistors  over  a limited 
voltage  range,  conventional  resistors  were  used.  The  contact 
size  was  approximately  0.03  in.  dia.  and  eight  heat  loss  resis- 
tors were  arranged  0.4  in.  apart  as  shown  in  Fig.  3*  The  values 
of  the  individual  resistors  were  determined  by  linear  approxima- 
tions to  the  electrical  equivalent  of  the  heat  loss  curve  over 
the  approximate  temperature  ranges  (derived  from  the  numerical 
solution).  The  appropriate  current-voltage  curve  was  obtained 
by  use  of  the  relation 

q _ 106  At 

i m A E 

where  A t = l44  F/volt,  as  measured. 

A E 


- 7 - 


These  resistor  values  check  very  closely  (approximately  6%) 
the  values  obtained  by  employing  the  "ef f ective"  surface  heat 
transfer  coefficient  of  Fig.  4 (Ref.  4),  and  using  the  relation 


1 . 
“hff 


4.  Discussion  of  Results 

Temperature  rise  data  on  the  two-dimensional  analog  model 
of  the  reinforced  concrete  T-beam  problem  were  obtained  at-  several 
locations  as  indicated  in  Fig.  3*.  For  direct  comparison  with  the 
numerical  results  given  in  ref.  [4],  all  data  were  taken  at  a 
thermal  time  of  2 hours  after  application  of  the  standard  tempera- 
ture-time fire  exposure  along  the  surface  designated  F.  Section 
A-A  is  sufficiently  removed  from  two-dimensional  effects  and  hence 
the  data  in  Fig.  4 represents  the  one-dimensional  temperature  gra- 
dient for  a horizontal  slab  whose  lower  surface  is  fire-exposed 
and  which  loses  heat  from  its  upper  ("unexposed")  surface  by 
Newtonian  cooling.  The  analog  results  are  within  7 %<>  based  on 
the  value  of  the  input  signal , of  the  corresponding  results  ob- 
tained by  numerical  solution  of  the  finite  difference  equations. 

In  Figures  5 and  6 are  shown  the  temperature  gradients  along 
Sections  B-B  and  C-C,  where  two-dimensional  heat  flow  patterns 
are  evident.  Here,  the  analog  and  numerical  values  do  not  differ 
by  more  than  11%. 

Higher  analog  values  along  part  of  the  B-B  axis  of  symmetry 
may  be  due  to  difficulties  associated  with  measurements  along  an 
edge  surface.  Although  this  problem  also  exists  along  the  other 
edge,  its  effect  is  partially  masked  by  the  Influence  of  the 
heat  loss  resistors  at  the  "unexposed"  surface.  Some  of  the  fac- 
tors involved  in  edge  measurements  include: 

(1)  the  conventional  edge  effect  for  a parallel  plate 
capacitor  which  results  in  increase  in  capacitance 
close  to  the  edge  [81. 

(2)  the  finite  dimension  of  the  probe  tip  and  its  prox- 
imity to  the  edge, 

and  (3)  the  imperfect  registration  of  the  grounding  plate, 
which  is  actually  recessed  approximately  0.03  in. 
from  the  edge. 
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To  avoid  the  edge  effect  at  the  axis  of  symmetry,  it  would 
be  preferable  to  utilize  the  complete  T-section.  This  would  also 
provide  twice  the  measuring  capabilities  (and  averaging-  advan- 
tages), but  would  probably  require  a larger  area  model  for  rea- 
sonable accuracy. 

In  most  experimental  work,  the  difference  or  error  between 
the  measured  and  true  potentials  is  expressed  as  a percent  of  the 
true  potential.  However,  In  analog-system  measurements,  this 
expression  is  not  suitable  since  the  true  potential  may  be  very 
small  or  even  zero.  A more  feasible  method  for  gaging  the  per- 
formance of  an  analog  system  is  to  express  the  difference  between 
the  measured  and  true  potentials  as  a percent  of  the  largest  po- 
tential difference  existing  in  the  field,  generally  the  voltage 
difference  between  boundary  electrodes  [2].  This  method  has  been 
employed  here  to  compare  the  analog  and  the  finite  difference 
results,  neither  of  which  can  be  considered  true  values. 

It  should  be  borne  in  mind  that  solutions  of  transient  heat 
flow  problems  by  finite  difference  computations  are  only  approx- 
imations and  involve  errors  which  depend  upon  the  space  and  time 
grids  chosen  and  the  associated  problems  of  instability  and  con- 
vergence. It  is  generally  difficult  to  estimate  the  magnitude 
of  such  errors.  Because  of  a lack  of  knowledge  of  the  accuracy 
of  the  numerical  computations,  analog  measurements  may  be  as 
much  a check  on  finite  difference  calculations  as  the  reverse. 

In  fire  endurance  tests,  the  performance  of  the  test  con- 
struction is  measured  by  the  time  period  during  which  it  with- 
stands the  standard  fire  exposure  without  failure  according  to 
prescribed  physical  and  thermal  criteria.  For  wall  or  floor 
constructions,  thermal  failure  generally  involves  a tempera- 
ture elevation  on  the  "unexposed"  surface  of  2^0  degrees  F. 
Structural  failure  of  stressed  structural  steel  members  usually 
occurs  at  a temperature  of  about  1000  F and  therefore  particular 
attention  was  paid  to  the  lower  portion  of  the  beam  where  rein- 
forcing members  would  be  most  subject  to  high  temperature. 

Studies  of  this  kind  are  useful  both  for  estimating  relative 
fire  endurance  times  and  as  an  aid  in  selecting  appropriate  steel 
and  concrete  combinations  for  particular  applications. 
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The  feasibility  of  using  a continuous-resistance,  continu- 
ous-capacitance analog  of  the  type  described  for  two-dimensional 
transient  heat  flow  analysis  has  been  demonstrated.  Considering 
the  size  limitations  of  the  electrical  model  and  the  variations 
in  its  electrical  properties,  the  comparison  with  the  numerical 
solutions,  at  least  for  the  locations  and  time  period'  measured, 
were  considered  reasonable.  The  exploratory  nature  of  the  ana- 
log model  described  in  this  report  leave  a number  of  questions 
regarding  experimental  technique,  measurement  error  and  overall 
accuracy  as  yet  unevaluated.  In  particular,  the  magnitude'*  of 
possible  errors  introduced  by  edge  effects,  by  contact  resis- 
tance and  by  capacitive  coupling  between  probes  and  lead  wires 
should  be  evaluated.  For  situations  requiring  increased  pre-ci- 
sion,  a more  thorough  exploration  should  be  made  of  the  dielec- 
tric properties  of  barium  titanate  (or  alternate  dielectrics), 
especially  with  regard  to  piezoelectric,  temperature,  hysteresis 
nonisotropic  and  dielectric  loss  effects. 
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Figure  l.\;  Schematic  diagram  for  measurement  of  the  time 
constant  of  a continuous-res  is tance,  continuous -capacitance 
analog. 


Figure  20  Mounting  details  of  analog  model 
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Figure  3.  Grid  pattern  for  analog  model 
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